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Abstract 

We present new results on the summation of nonpower-series expansions 
in QCD Analytic Perturbation Theory (APT) in the one-loop approx- 
imation. We show how to generalize the approach suggested by one of 
us earlier to the cases of APT and Fractional APT (FAPT) with heavy- 
quark thresholds. Using this approach we analyze the Higgs boson decay 
H — > bb. We produce estimations of the higher-order corrections im- 
portance in (F)APT and present a very transparent interpretation of the 
resummation results in terms of Aqcd shifts. 



1 Introduction 

This is a written version of the talk given at the Memorial Igor Solovtsov Semi- 
nar, held on January 17-18, 2008, at the Bogoliubov Laboratory of Theoretical 
Physics (BLTPh), JINR (Dubna). For this reason first we want to recollect our 
contacts with Igor Solovtsov in and outside the lab. 

We met first time at BLTPh in the mid 90s when he conducted a series 
of seminars on the variational approach to quantum field theory, in particular 
to QCD. He was those times very enthusiastic about this approach which later 
on guided him to Analytic Perturbation Theory. Only after some time other 
people started to share his enthusiasm and understand what he knew at the very 
beginning. We met ones at the Institute for Theoretical Physics of Heidelberg 
University (presumably, in 1997) and participated in the seminar on Analytic 
Perturbation Theory he presented there — it was a kind of enjoinment, so clearly 
he presented the matter. We remember very well how he was preparing his 
doctorate thesis — calmly, with some kind of dignity, and remaining to be very 
open for discussions. Later on, we met also during some conferences. The one we 
remember more clearly, is the Research Workshop on Calculations for Modern 
and Future Colliders, held in Dubna in July 2000. The reason for this memory 
is especial: during the boat-trip along the Volga river it appeared that Igor was 
a very good singer — he had a good voice, an absolute pitch and a very nice 
song repertory. One of us, A.B., met Igor last time quite recently, during the 
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XIV International Seminar "Nonlinear Phenomena in Complex Systems" (May 
22-25, 2007, Minsk) and ones again enjoyed his nice voice. To resume, Igor was 
a talented physicist, a nice man, and a good friend. Let God bless his memory 
and Igor be remembered forever. 

Now we continue with his 'child' — APT. Since the original works of Jones 
and Solovtsov [1,2] and Shirkov and Solovtsov [3] appeared in 1995-97, the 
analytic approach to QCD perturbation theory has progressed considerably. The 
main object of this approach is the spectral density that provides the means to 
define the analytic running coupling in the Euclidean region using a dispersion 
relation. The same spectral density is used to define also the running coupling 
in Minkowski space, by employing the dispersion relation for the Adler function 
[4,5]. The analytic approach was extended beyond the one-loop level [5,6] and 
special analytic and numerical tools were elaborated [7-9]. APT appeared to 
be well-suited to deal with more than one (large) momentum scale [10, 11]. All 
these efforts resulted in a systematic approach, named Analytic Perturbation 
Theory (APT), recently reviewed in [12]. 

In order to treat hadronic amplitudes in the same manner one needs to 
use not only analytization of the running coupling, but also more complicated 
objects, like integrated couplings J^dxj^dy a s (Q 2 xy) f(x)f(y), as was firstly 
indicated in [13, 14]. The authors of this paper proposed the principle of ana- 
lytization as a whole to treat such objects and new techniques was developed 
to produce needed analytic images. The same principle applies also to other 
hadronic amplitudes in the renormalization-group based QCD factorization ap- 
proach, where complicated objects, like fractional powers of running coupling 
and its products with logarithms, ln(Q 2 //Li 2 ), naturally arise. A version of the 
"analytization" approach, which takes into account all terms that contribute to 
the spectral density, was developed in [15-17] and named Fractional Analytic 
Perturbation Theory (FAPT). 

In this paper, we present new results on the resummation of non-power- 
series expansions in APT in the one-loop approximation. We show how to 
generalize the approach suggested by one of us earlier [18] to the cases of: (i) 
global (with taking into account heavy-quark thresholds) APT, and (ii) global 
FAPT. Using this approach we analyze the Higgs boson decay into a 66-pair of 
quarks. We produce estimations of how important higher-order corrections can 
be and present very transparent interpretation of the resummation results in 
terms of Aqcd shifts. 

2 Fixed-TV/ and global APT 

First, we need to explain our notation. In order to have a direct connection to 
our previous papers [15-17] and to simplify the main formulae in the sections, 
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where we speak about fixed-order (F)APT with fixed value of active flavors 
Nf, we use here the normalized coupling [18] a(Q 2 )=(3f a s (Q 2 ; Nf) and analytic 
images are constructed for powers of this object: 

A n {Q 2 ) = AeK] = J^J^Lda; (la) 
2l„( S ) = AmK] = r^-da (lb) 



with p n (a) = Im [a n (-a)} and f = b (N f )/(4ir) = (11 - 2N f /S)/(4ir) with 
bo(Nf) being the first coefficient in the QCD (3 function. 

When we discuss a global version of (F)APT, where Q 2 (or s) varies in the 
whole domain [0, oo) and Nf effectively becomes dependent on Q 2 (or s), we use 
the standard running coupling a s (Q 2 ), and analytic powers are defined for the 
corresponding powers of this coupling 

fOO glob/' \ 

■Ai oh (Q 2 ) = A E [a£] = / ^4 da; (lc) 

Jo & -r v 

/•OO glob/ \ 

2lf b ( S ) = A M K] = / da . (Id) 

In order to distinguish these two approaches, we introduce the upper index glob . 
Here we need also analogs of fixed- Nf quantities with standard normalizations: 

MQ 2 ;N f ) = ^l; % n (s;N f ) = ^, (le) 

which correspond to analytic couplings A n and 2l n in the Shirkov-Solovtsov 
notation [12]. These couplings also have integral representations of the type 
(la)-(lb) with the spectral densities p n (a;Nf) = p n (a)//3j. 

The spectral density p^ lob (o") is simply related with the spectral density 
p n (a;N f ) [3,7,19]: 

6 

pf » = 6 [a < m 2 ] p n (a; 3) + £ 9 [m 2 < a < m} +1 ] p n (a; Nf) , (If) 

/=4 

where, in order to make our formulae more compact, we denoted m& = m c , 
= mj, uiq = nit, and ra-j = +oo. 

We need also to speak about effective coupling as a function not of Q 2 or s, 
but rather of logarithms L = ln(Q 2 /A 2 ) or L = ln(s/A 2 ). In these cases, we use 
the same notation for coupling but with an argument placed in squares: a"[L], 
*4j,[L], and 21^ [L]. Then, in the one-loop approximation we have 

*[L] = l "M = Z2Ti*' * [L > N ' ] = p f [V + **\ - (2) 
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2.1 Series summation in the one-loop APT with fixed Nj 

Let us consider different sorts of perturbative series expansions of a typical phys- 
ical quantity, like the Adler function, D[L], and write at the one-loop level [17], 




do + di ^2 d~v 



n=l 



a n [L] 

A n [L] 



(3) 



where L = In (Q 2 /A 2 ) applies to the Euclidean {D[L], V[L], a n [L] and A n [L\) 
and L = In (s/A 2 ) to the Minkowski (TZ[L] and 2l n [L]) regions, and d n = d n /d\. 
It is useful to introduce a generating function P(t) for the series expansion whose 
specific form depends on a particular quantity in question: 



d n = P(t) t n ^dt with / P(t) dt = l. 
Jo Jo 

To shorten our formulae, we use a shorthand notation 

POO 

{{f(Jt)))p= / f(t)p(t)dt, 

Jo 



(4) 



(5) 



where the brackets ((. . .))p denote the average over t with the weight P(t). 

The different couplings a, Ai, and 21 1 satisfy the same one- loop renormaliza- 
tion-group equation that can be rewritten as a one-loop recurrence relation [20- 
22] 



7 n+l 



[L] 
A n+1 [L] 




(6) 



Substituting Eqs. (6), (4) into the perturbative series expansion, Eq. (3), one 
obtains 



f D[L] 
V[L] 

( n[L] J 



> = do + dxY^ 



«(-tr»p d n 



n=0 



nl 



dL n 



a[L] 
Ai[L] 
ML] J 



= do + di 



((a[L-t])) P 
({Ai[L-t])) P 

«ai[L-t])) P 



(7) 



as it was shown by one of us (S.M.) in [18]. As long as we have not proved that 
the order of summation and integration can be interchanged, this representation 
has only a formal meaning. Note, however, that the integration over the Taylor 
expansion of the term a[L — t] in the integrand reproduces the initial series for 
any partial sum. It is not surprising that the integrand in the standard case (first 
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line in Eq. (7)) faces a pole singularity and is, therefore, ill-defined. On the other 
hand, for the last two cases the integral has a rigorous meaning, ensured by the 
finiteness of the couplings, i.e., Ai[t] , 2li[i] < 1. Since any coefficient d n is the 
moment of P(t), this function should fall off faster than any power, e.g., like an 
exponential or faster. Therefore, the APT expressions on the RHS of Eq. (7) 
exist and are proportional to A\[L — i{L)] or 2li[L — i(L)], where i{L) can be 
treated for each L as an average value of t associated with this quantity. 

Provided the generating function P{t) is known, one can compute the aver- 
age ((Ai[L — t]))p (or — t]))p) in Eq. (7) explicitly and obtain all-order 
estimates for the expanded quantity under consideration. Unfortunately, only 
the first few d n coefficients are known for most relevant processes. However, 
the asymptotic tail calculated first by Lipatov et al. [23], see for review [24], 
allows us to construct a model of the generating function P(t) with reasonable 
accuracy: 



d n ~ r(n + 1) n^c n 



1 + 



P(t) ~ (Vc) 7+1 e-*/ c , (8) 



where 7 < 1 and c are numerical coefficients. 

Now we illustrate this statement by the example of perturbative coefficients 
d n for the Adler function 



D S [L] = d. + d^dJ^X 

n=l ^ ' 



(9) 



for scalar correlator, cf. [17,25], see the first row in Table 1. To simulate the 
n-dependence of the coefficients d n , in accordance with recipe (8), we use the 
model: 

g = ^, r(» + i) ++ w») [ (Wa) 

which is produced by the generating function 

In the second row of Table 1 we show the results obtained by fitting two known 
coefficients d2 and ds using the model (10) 1 . We see that it gives a very good pre- 
diction for $1 = 662, which appears to be quite close to the value 620 calculated 
by Chetyrkin et al. in [25] a year ago. If we use instead the fitting procedure 
with taking into account the fourth coefficient d^ then the parameters c and (3 of 
our model (10) change slightly: {c = 2.5, (3 = -0.48} -> {c = 2.4, (3 = -0.52}. 
This shows the reasonability of our modelling. 



1 Note that is automatically equal to unity. 
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Model 


di d 2 d 3 di d 5 


pQCD results with TV/ = 5 


1 7.42 62.3 620 


Model (10) with c = 2.5, f5 = -0.48 


1 7.42 62.3 662 8615 


Model (10) with c = 2.4, (3 = -0.52 


1 7.50 61.1 625 7826 


"NNA" prediction of [26] 


1 3.87 21.7 122 1200 



Table 1: Coefficients d n for Higgs boson decay series. 



Predictions for d n values in the so-called "Naive Non-Abelianization" (NNA) 
approximation [26] are presented in the fourth line of Table 1 and they appear 
to be significantly smaller than the exact values as well as our model estimates. 



2.2 Series summation in the global APT 

Consider the situation with resummation in the global APT when we take into 
account heavy-quark thresholds. We consider here the following series: 

oo 

TZ^ h [L] = d Q + d 1 ^J n 2lf b [L]. (11) 

n=l 

Note that due to different normalizations of 2l„[-L] and 2l^ lob [L] the coefficients 
d n in Eqs. (3) and (11) are different. In Appendix A we consider in detail the 
case with only one threshold, corresponding to the transition Nf = 3 — > Nf = 4, 
and provide the result of summation in the Minkowski region in Eqs. (A6)-(A7). 
With taking into account all other thresholds the final result reads 



6 r t i 

TZ^ h [L]=d + dx J](9(L / <L<L /+1 )«2ii L+X f -—;f » P 

/=3 Pf 

5 6 

+ d 1 Y,HLf<L<L f+1 ) Y, «A fc ai[t]» P , (12) 

/=3 k=f+l 



where Aj = ln(Ag/Aj) describes the shift of the logarithmic argument due to 
the change of the QCD scale parameter (A3 — > A/), Lf = h^m^/Ag), L3 = — 00 
and Lj = +00. The second term in Eq. (12) looks like a natural generalization 
of the fixed- Nf formula (7) with taking into account different QCD scales Af for 
each fixed- Nf region, while the last term appears due to continuity of 7£ glob [I/] 
at heavy-quark thresholds. Using a toy model Ptoy(*) = S(t — 2) we estimate 

6 

a relative contribution of the J2 ((Afc2lx[i]))p-term: it is of the order of 3% 

k=f+i 
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for the / = 4 threshold, 0.5% for the / = 5 threshold, and 0.1% for the / = 6 
threshold. 

In the Euclidean region the result of summation, as shown in [27], is more 
complicated: 

f=3 JL f ~ l ~ 
6 

+ d 1 J2((^f[L,t]))p, (13) 
where corrections ((Af[L,t]))p to the naive expectation formula are defined as 

- /' ftl£/+ y-'-^-:'V'^. (i4) 

Jo I3f-i [l + e L - L l-"=^i-<] 

and now, in contrast to the Minkowski case, they explicitly depend on L. To 
estimate the effect of the thresholds in the Euclidean domain, we compare the 
numeric value of (13) for our toy model ftoy(i) = 5(t — 2) with the numeric value 
of the 'naive' summation formula 

The relative values of the differences between (13) and (15) are of the order of 
2.5% at L ~ 5 and are about 0.5% for L > 12. 

3 Fixed-vV/ and global FAPT 

Sometimes perturbative series starts not with the unity term like do in Eq. 
(3) but rather with doa u [L] with power v being a fractional number, see, for 
example, [17]. Then one needs to use FAPT in order to obtain an analytic 
image of the series. For this reason we need to generalize the APT summation 
method discussed in section 2 in order to apply it in FAPT case. 

3.1 Resummation of expansions in the fixed- Nf FAPT 

First, we consider the case of fixed Nf. We start with a power series in the 
standard perturbative QCD 

oo 

D V [L\ = do a u \L\ + d 1 J2 a n+v [L] , (16) 
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with v > that appears due to renormgroup summation. This type of series 
generates the following analytic image in the Minkowski region: 



n=0 



K v [L] = d % v [L] +d!^2 dn+l %n+l+u [L] . 

Due to the recurrence relation 
we have 



i» 



? \ n 

4 **m 



T{n + v) V dL / 
K U [L] = d Vl u [L}+d 1 ((X(t;l + v))) P , 



(17) 



(18) 



(19) 



where 



x{t-i+v) = ^ ! ,:, r/ , i " l ', (l \" ] * / win. .*•,-/ — . (20) 



n=0 



T(n+ 1 + v) 



dL 



We have a nice representation for the ratio of T-functions (see formula 5.2.7.20 
in [28]) 



r(n+l)I> + l) 



T(n + u + l) 

which gives us a possibility to process further with our series (20) 
i-x t ) n r(l + v) 



(21) 



E 

n=0 



F(n + 1+u) 



Now we recall that xt is the operator xt = t-d/dL. That means that the operator 
e zx t w h en acting upon a function -A[L] just shifts its argument: e ZXt A[L] = 
A[L + iz], and we have 



X(t; 1 + v) = J % 1+v [L + t (V^ - 1 



du . 



(22) 



Substituting Eq. (22) in Eq. (19), we obtain as a result: 



K U [L] = d 2UL] + di du{{VL 1+v L-t 

Jo 1 
= d 2l,[L] + d 1 ((2t 1+ ,[L-t])) Pi/ , 



(23) 
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where 

$„(z) = vx v ~ x — ► 6(x). (24b) 

Note here that lim P v (t) = P(t). 

A similar formula can be obtained to sum up perturbative series in the Eu- 
clidean domain: 

V„[L] = (ioAW + d 1 ((A+,[i-t]))p,, (25) 

We see that the main difference between FAPT and APT with respect to 
a series summation technique is that in FAPT one needs to modify the initial 
generating function P(t) into a new one, P v (t). It is interesting to note that for 
our model 

*• = "-^y = (26) 

the integration for P u (t) in Eq. (24a) can be carried out exactly to produce rather 
a complicated expression containing the regularized confluent hypergeometric 
function. For integer values of 5 = m > formula is simplified to 



v 
0, m 



P„(t;m) = T(v+l)G 2 {° 2 
where G^'j^l...) is the Meijer G-function defined as 



(27) 



n 2 <° 

^1,2 



h, b 2 J 2-Ki J c T(a + s) 



(28) 



where the contour C of integration is set up so that the poles of T{b\ + s) and 
r(&2 + s ) ne in the same region (internal or external) with respect to the contour 
C. 



3.2 Resummation of expansions in the global FAPT 

To generalize formulae obtained in the previous subsection to the case of the 
global FAPT, we follow the lines of subsection 2.2 and take into account the 
new recurrence relation (18) for couplings in FAPT. We start with the analytic 
image of perturbative series in the Minkowski region 

oo 

Kt° b [L] = do^m+d^dn+i^tti+AL} (29) 

?1=0 
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and proceeding carefully obtain the result 2 : 



TZf b [L] = d K loh [L] +d l Y,0(L f <L< L f+1 ) ((% 1+v [l+ X f - — ; /] » P „ 

f=3 ' 1 f 

6 6 

+ d 1 Y J 0(Lf<L<L f+l ) Y, «A fc 2ti + „[t]»p„. (30) 

f=3 k=f+l 

Formulae for the Euclidean region appear to be more complicated and we refer 
the interested reader to the review [27], where this case is analyzed and all 
formulas are derived. 



4 Application to the Higgs boson decay 

In order to apply the developed techniques to the estimation of real processes, 
we need first to model the generating function P(t) of the corresponding per- 
turbative coefficients. As we know from the discussion in section 2, in the case 
of the Higgs boson decay H° — ► bb the model (10) is of reasonable accuracy. 

First, we recall what is the object of interest in the case of the Higgs boson 
decay to a bb pair. Here we have for the decay width 



T(H bb) 



G F 



M H R S (M 2 H ), 



(31) 



where Rs(M^) is just m^(M^) Rs(M]j). In the one-loop FAPT this generates 
the following non-power expansion 3 : 



pFAPT 



[L] = 3m 



(i) 



n>l 



(32) 



where is the renormalization-group invariant of the one- loop mj(/i 2 ) evolu- 
tion 



m 2 b (Q 2 )=m 2 {1) a^(Q 2 ) 



(33) 



with vq = 27o/6o(5) = 1-04 and 70 is the quark-mass anomalous dimension. 

We take the model (10) and apply Eq. (30) to estimate how good is FAPT 
in approximating the whole sum Rl APT [L] in the range L G [11,13.8] which 

corresponds to the range M H £ [60, 170] GeV 2 with = 172 MeV and 



2 With the same notation L3 = —00 and 1/7 = +00 as in Eqs. (12) and (13). 

3 Appearance of denominators it™ in association with the coefficients d n is a consequence of 



dn normalization, see Eq. (9). 
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Fig. 1: The relative errors A.n[L], N = 2,3 and 4, of the truncated FAPT, Eq. 
(35), as compared to the exact summation result, Eq. (34). 



2lf lob (m|) = 0.120. In this range we have L 5 < L < L 6 , so that Eq. (30) 
transforms into 



pFAPTrri j 



3m 



(i) 



7T 



L+X 5 — ^-;5 



»P, n (34) 



with P UQ (t) defined via Eqs. (10) and (24a) with the parameters c = 2.4, (3 
-0.52, and v = 1.04. 

Now we analyze the accuracy of the truncated FAPT expressions 



i? FAPT [L;iV] = 3 m?! 



(i) 



JV 



n=l 



(35) 



and compare them with the total sum i?g APT [L] in Eq. (34) using relative errors 



nFAPT 



[L;N] 



pFAPT 



[L] 



(36) 



In Fig. 1, we show these errors for N = 2, N = 3, and N = 4 in the analyzed 
range of L G [11,13.8]. We see that already i?| APT [L;2] gives accuracy of the 
order of 2.5%, whereas J R| APT [L; 3] of the order of 1%. That means that there is 
no need to calculate further corrections: at the level of accuracy of 1% it is quite 
enough to take into account only coefficients up to d%. This conclusion is stable 
with respect to the variation of parameters of the model Pn(t). These estimates 
demonstrate also a good convergence of the considered series in FAPT. 

It is also interesting to estimate how large is the effective shift AL which 
models the exact result (34) by the simple formula 



R 



FAPT 



L) 



o - 2 

3m (1) 



Of glob 
+ 



L- AL 



(37) 
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Numerical estimation shows that AL ~ 2.16. We can transform this result into 
a finite shift of Aq£= 3 = 0.172 MeV: 

A QCD 3 - A e ^= 3 = 2.94 A^ 3 ~ 505 MeV . (38) 

If one estimate this shift by using only the cfe-term, then the result is 10% 
smaller: A cff /_ ~ 450 MeV. 

5 Conclusions 

In this paper, using recurrence relations of the one-loop (F)APT, see Eqs. (6) 
and (18), we construct all needed rules to sum up nonpower series expansions 
in the fixed- Nf and global (F)APT. We show how these results can be used to 
estimate the rate of convergency of nonpower series expansion. In the case of 
the series appearing in the QCD description of the Higgs boson decay into a 
56-pair, this convergence is estimated to be moderately fast: in order to have 
1%-level of accuracy, it is sufficient to sum up the first four terms (do, d±, d,2, 
and ds). In this case, such summation is equivalent in the considered range of 
L e [11, 13.8] to the constant shift of A^ 3 from 172 MeV to 505 MeV. 
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A Single quark threshold in the Minkowski APT 

We consider here only one heavy-quark threshold corresponding to the transition 
Nf = 3 — > Nf = 4. In this case, the global one-loop spectral density pfl° h (s) = 
Pn° b {L(s)] (with L(s) = In (s/A|)) is expressed in terms of fixed-flavor spectral 
densities with 3 and 4 flavors, p n [L; 3] and p n [L + A4; 4], see Eq. (If): 

pC h [L] = 9(L<L 4 ) p n [L;3] + 9(L 4 <L) p n [L + X 4 ;A] , (Al) 
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with A4 = In (A|/A|) describing a shift of the logarithmic argument due to the 
change of the QCD scale parameter (A3 — ► A4), L 4 = In (m 4 /A|), and p n [L; Nf] 
defined in Eq. (1). In the one-loop approximation we have a very useful property 
of the spectral densities p n [L] and p n [L; Nf]: 



Pn+l[L] 

p n+1 [L;N f ] 



- ( 4r) Pn[L] 



n 



d_ 

dL J r " L ~ J r(n + l) V dl 
d 



n(3 f 



dL 



Pn\L;Nf] 



n+l 



1 \ n 

Tl) 



— ) Pl [L]; (A2a) 
(A2b) 
(A2c) 



which is valid for n > (except only Eq. (A2b)) and allows us immediately 
to rewrite Eq. (Id) in a more explicit form, where all n-dependencies appear 
explicitly: 



An+l[ 1 ~ F(n + 1) 1/3? 



+ 



(n+l) 

1 

+ M 

9(L> U) 1 
r(ra+l) 7% 



.±) n ^[L;3]- (-A-YmU-A 
d 



dL^ 

d^ n 



dL± 

2ti[L 4 + A 4 ;4]} 



-— J 2li[L + A 4 ;4] 



(A3) 



Thus the general structure of the n-dependence, which we have in Eq. (A3), is 
x n /r(n + l) with x = [l//3/](— d/dL) and we know how to sum up this structure, 
see Eq. (7): 

00 

n *<*> [L] ^ do + J2 d n+1 [L] ^do + d^ Oi [L] 6 f ,i [L + Xf] ; (A4) 



n=0 



/3?r(n + i) V^l y 



=0^7 



j\ n 

-£) a 1 [L;JV>] = «Si 1 



L — — : N 



)>p. (A5) 



Collecting all types of different ^-structures in Eq. (A3) and inserting them in 
Eq. (A5) we arrive at the answer 



K gloh [L] =d + d 1 6(L< L 4 ) <<2ti 
+ d 1 6(L> L 4 )((2li 



P3 - 



+ A4&1 [t]»p 



L + A 4 --f ;4 )) P , 

P4 



where we denoted (with A3 = 0) 

Af* n [t]=* n [Lf + \f-±;N f 



21, 



0/- 



/-1 



(A6) 



(A7) 
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B Numerical details 

We use the following heavy-quark masses: m c = 1.2 GeV, m& = 4.3 GeV, and 
m t = 175 GeV that generate the following set of QCD scales and logarithms 
Xf = In Ag/Aj in the one-loop approximation: 



A, 



172 MeV, A 4 = 147 MeV, A 5 = 111 MeV, A 6 = 55 MeV; (Bl) 



A 4 = 0.311, A 5 = 0.898, A 6 = 2.30, \ z = In 



12.55. 



(B2) 
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